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Abstract 

The explicit violation of the general gauge invariance/relativity is adopted as the 
origin of dark matter and dark energy of the gravitational nature. The violation of 
the local scale invariance alone, with the residual unimodular one, is considered. Be- 
sides the four- volume preserving deformation mode - the transverse-tensor graviton 

- the metric comprises a compression mode - the scalar graviton, or the systolon. 
A unimodular invariant and general covariant metric theory of the bimode/scalar- 
tensor gravity is consistently worked out. To reduce the primordial ambiguity of the 
theory a dynamical global symmetry is imposed, with its subsequent spontaneous 
breaking revealed. The static spherically symmetric case in the empty, but possibly 
for the origin, space is studied. A three-parameter solution describing a new static 
space structure - the dark lacuna - is constructed. It enjoys the property of gravi- 
tational confinement, with the logarithmic potential of gravitational attraction at 
the periphery, and results in the asymptoticaally flat rotation curves. Comprising a 
super-massive dark fracture (a scalar-modified black hole) at the origin surrounded 
by a cored dark halo, the dark lacunas are proposed as a prototype model of galax- 
ies, implying an ultimate account for the distributed non-gravitational matter and 
a putative asphericity or rotation. 

PACS numbers: 04.50.Kd Modified theories of gravity - 95.35.+d Dark matter 

- 95.36.-I-X Dark energy 



1 Introduction and motivations 

The General Relativity (GR) is well-known to be the viable theory of gravity perfectly 
consistent (modulo some additional assumptions) with all the available observations. 
Nevertheless, it hardly is an ultimate theory being rather an effective field theory pro- 
duced by a more fundamental one. The latter well may enjoy a wider "low-energy" 
remnant, with GR being just a principle part of it. The basic requirements to such 
a GR extension may be that it should safely retain all the well-established theoretical 
and observational properties of GR, but at the same time should encompass some new 
phenomena beyond GR (or v. v.). Consider these items in more detail. 



Gauge invariance: unimodular vs. general The essence of GR may be expressed 
by saying that it is a gauge theory of a massless tensor field/graviton. As the respective 
gauge group it is conventionally taken the group of the general diffeomorphisms. For 
consistency with GR, an extension to the latter should conceivably be based on a gauge 
principle, too. Following the original description of the irreducible (one-dimensional) 
unitary representations of the (non-compact) Poincare group due to Wigner [1], it was 
later found by van der Bij, van Dam and Ng [2] that the necessary and sufficient gauge 
group admitting a massless tensor field is the group of transverse diffeomorphisms, with 
the group of general diffeomorphisms being thus excessive to this purpose. For this 
reason, a theory of gravity based on such a minimally violated gauge invariance might 
well be the most natural candidate to supersede GR. To retain general covariance the 
transverse diffeomorphisms are to be substituted though by the unimodular ones (see 



To clarify the structure of such a GR extension decompose the group of the general 
diffeomorphisms into the commuting subgroups of the transverse diffeomorphisms and 
that of the local scale transformations. By construction, the determinant of metric, 
(yf, changes only with respect to the latter subgroup. Under restriction by the transverse 
subgroup, g is a scalar and can be treated as an independent field variable to be implement 
in the gravity Lagrangian threefold. 

(z) General invariant Lagrangian First, one may retain g in the gravity Lagrangian 
implicitly through metric, but made it unphysical by means of the local scale transfor- 
mations as in GR (or in its modifications preserving general gauge invariance). As a 
physical degree of freedom in metric there is left just the (massless) tensor graviton. 

(a) Restricted Lagrangian Second, one may restrict the Lagrangian by fixing g a 
priori (as it was originally proposed in other context by Anderson and Finkelstein [3]), 
so that the respective degree of freedom in metric is absent, leaving only the massless 
tensor graviton (see also, e.g., [2], |1]-[I7]). Being a close counterpart of GR such a 
theory of tensor gravity is conventionally referred to as the Unimodular Relativity (UR). 
A possible advantage of UR compared to GR is that the cosmological term is no more a 
Lagrangian parameter as in GR itself, but appears as an integration constant. This may 
help in solving the so-called naturalness problem of cosmological term. 

[Hi) Extended Lagrangian And v. v., one may extend the gravity Lagrangian in the 
wake of Buchmiiller, Dragon and Kreuzer [IH]^[2n] by adding to it the terms explicitly 
containing the derivatives of g. Due to violation of the general gauge invariance, the 
degree of freedom corresponding to g can not be eliminated any more and gets physical. 
Thus, in addition to the massless tensor graviton the metric comprises a (light) gravita- 
tionally interacting scalar particle. Studying such a scalar-tensor extension to GR is of 
considerable interest by itself from the theoretical point of view. 

Unimodular invariance and dark matter However, the main motivation to adhere 
to the theory of the latter type may lie in the problems of the so-called dark matter 



Universe, present perhaps the greatest challenge facing the modern physics. Basically, 
while DE should be "spilled" all over the Universe on the cosmological scale, DM should 
cluster on the galactic and galaxy cluster scales. At that, DM proves to be one of the 
most important building-blocks of the Universe. Its energy yield largely exceeds that of 

^For a review on these topics see, e.g., [21]. 



Sec. 2.1). 




substances, being crucial for the structure of the 
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luminous matter. While DE might already be anticipated as given by the well-known 
cosmological term, DM appears, surprisingly enough, absolutely ad hoc from the modern 
theory point of view. Although the evidence for DM is overwhelming, its nature remains 
still obscure. 

The crucial points concerning DM are, first, that it should interact with the lumi- 
nous matter very faintly and, second, its clusterization in gravity field should be weaker 
compared to the luminous matter. Such an elusive kind of matter may well have an un- 
conventional nature, e.g., the gravitational one. With this in mind, one can put forward 
the hypothesis that as an origin of such a DM there may serve the explicit violation of 
the general gauge invariance/relativity [22] |^ The reason is that the ensuing in this case 
extra terms in gravity equations may be treated as an effective energy-momentum tensor 
of some additional physical degree (s) of freedom. Interacting only gravitationally, but 
otherwise than the non-gravitational matter, this new substance may well suit as DM. 
Thus due to violation of general relativity the metric field itself may store a kind of the 
gravitational DM|^ 

As a paradigm one may consider the minimal violation of general relativity to the 
unimodular one [22l [23]. In this case, on the one hand, the residual invariance suffices 
to justify the masslessness of tensor graviton. On the other hand, due to violation of 
the local scale invariance alone, the metric comprises no more than one extra physical 
component besides the massless tensor graviton. Originating from metric and interacting 
only gravitationally, the respective particle may naturally be associated with the scalar 
graviton. To have sense in the arbitrary observer's coordinates the theory should be cast 
to the general covariant form by introducing a non-dynamical scalar density of the same 
weight as y/—g- By this token, it was found appropriate to treat the scalar-graviton 
field as an independent variable, embedding within it (virtually) the unknown scalar 
density. In the static spherically symmetric case this made it possible to find an exact 
solution to the extended gravity equations in the empty, except for the origin, space 
Such a solution describes a two-parameter generalization of the one-parameter black 
holes (BHs) and is valid also in GR with a (free massless) scalar field. Beyond GR, an 
approximate, regular at the origin one-parameter solution in the empty space was then 
found analytically [25] and refined numerically [26j. This solution presents a halo- type 
static space structure built entirely of scalar gravitons. Such a structure was shown 
to possess a soft-core energy density profile, qualitatively compatible with that for the 
galaxy DM halos. This makes it urgent to investigate the more general solutions to the 
aforementioned equations, as well as consequences thereof. 



Content In this article the metric theory of the scalar-tensor gravity built on principles 
of the unimodular relativity and general covariance is consistently worked out. In Sec- 



tion 2 the theoretical background is developed. In Sec. 2.1 the explicit violation of general 



relativity is discussed in toto. The violation of the local scale invariance alone, with the 



residual unimodular invariance, is then investigated in more detail. In Sec. 2.2 the ef- 
fective field theory of metric, enjoying the residual unimodular invariance as well as the 
general covariance, is considered. To terminate a priori allowed effective Lagrangians, 
a dynamical global symmetry at the classical level is imposed. The ensuing classical 



^The term "relativity /invariance violation" is more appropriate (see Sec. 2.1) than the "covariance 
violation" used previously in [52]. 

This does not exclude, of course, a fraction of a more conventional particle DM, so that the direct 
searches for DM are by no means meaningless. 
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equations are then written down, and the spontaneous breaking of the global symmetry 
is displayed. Section 3 deals with the ensuing classical equations and their solutions in 
the static spherically symmetric case. In Sec. 3.1 the gravity equations in the empty, 
but possibly for the origin, space are presented in a specific gauge. In Sec. 3]2] an exact 
two-parameter solution of the BH-type, valid also in GR with a (free massless) scalar 
field, is rebuilt in the most transparent fashion, and then compared in several gauges. In 
Sec. 3^ an approximate, regular at the origin one-parameter halo-type solution, missing 
in GR with a free scalar field, is exposed. In Sec. 3.4 an irregular at the origin three- 
parameter solution interpolating between the two preceding extreme cases is constructed. 
The property of the gravitational confinement for the respective extended static space 
structures is found. Section 4 is devoted to interpretation and application of the found 
solutions/structures in the context of galaxy DM. In Sec. 4.1 the energy content of the 
structures is revealed. In Sec. 4.2 the asymptotic flat rotation curves (RCs) ensuing due 
to these structures are exposed, with their DM interpretation presented. Finally, the 
relevance of the structures to galaxies is discussed. In Conclusion the state of affairs of 
the theory and its future prospects are outlined. 



2 Theoretical frameworks 
2.1 General relativity violation 

Relativity /invariance vs. covariance To go beyond GR let us reflne the paradigm of 
general relativity and/or general covariance adopted in GR. Let I = J Cd^x be the action 
of a fleld theory, with C being the density of its effective Lagrangian. To be precise, 
under the (classical) dynamics there will be understood the principle of the minimal 
action, with the emerging classical equations. In this respect, C may depend on two 
kinds of fields: the dynamical / relative and non- dynamical / absolute ones, generically, (^'s 
and (/?*'s. The non-dynamical fields are given a priori. The dynamical ones are those, 
classical equations for which are obtained by extremizing /, under frozen </9^,'s. The <^'s 
include metric g^i, (or its restricted part), matter fields and, optionally, the undetermined 
Lagrange multiplier which is a kind of the dynamical scalar field. The classical solutions 
for ip^s should depend on y9*'s as on the external functional parameters. An arbitrary 
effective field theory, in contrast to GR, may thus be said to be the theory of "restricted 
relativity", with y^'s implemented in the relative fashion and (y9^,'s in the absolute one. 

Accordingly, there are envisaged two types of space-time properties of the field theory. 

{i) Covariance This is a kinematic property describing a maximal kinematically 
allowed set G of the simultaneous diffeomorphisms of y^'s and v?*'s, under which I remains 
invariant. For consistency, (p^^s should transform under G through themselves, whereas 
the transformed yj's may depend on ip^^s as well. To have sense in the arbitrary observer's 
coordinates (chosen a priori) C should be a general covariant scalar density due to a 
required number of (^^,'s (if any). In what follows, G will be supposed to be the group of 
general diffeomorphisms. 

[a) Relativity This is a dynamic property describing a maximal subset if C G of the 
diffeomorphisms of y^'s alone (with y9*'s transforming as scalars), under which / remains 
invariant. At that, in the process of variation any 6ip within the configuration space are 
admitted, whereas (y9^,'s are to be frozen {6ip^: = 0). By the very nature, it is H which 
serves as a space-time gauge group for the theory. So, as a synonym of relativity there 
will also be used the term (gauge) invariance. 
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In GR the covariance and relativity/invariance coincide, both being the general ones. 
Beyond GR one should distinguish them. In particular, the non-trivial quotient G/H 
may result in appearance of the extra physical degree(s) of freedom. In what follows we 
consider the unimodular invariant but general covariant metric theory of gravity, with a 
dynamical scalar-graviton field in metric and a non-dynamical scalar-density field. This 
is the maximal allowed explicit violation of general relativity/invariance consistent with 
the masslessness of transverse-tensor graviton. 



Unimodularity and bimodality To be precise, under the theory of gravity at the 
energies lower than the Plank mass there will be understood the effective field theory 
of metric. Ultimately, the latter theory is to be considered as a "low-energy" remnant 
of a more fundamental theory. This remnant may basically be characterized by the two 
ingredients: a residual invariance group H C G and a set of the low-energy fields. As a 
paradigm, let H be the maximal subgroup of G given by the unimodular group U. Let 
Qfj^i, be a dynamical metric field, with det{g^i,) = g, and let ^/—g^ be a non-dynamical 
scalar density of the same weight as y/—g- Related to DM and associated with a non- 
dynamical measure, \/—g* will be, in the wake of [3], referred to as the (dark) modulus. 
A priori there are envisaged arbitrary (including singular) moduli. The dark moduli are 
to be considered as a kind of a new substance, the nature of which in the framework 
of the effective field theory is unspecified and should ultimately be revealed by a more 
fundamental theory. 

Without loss of generality, the Lagrangian density C for a unimodular invariant ef- 
fective field theory of metric may be expressed through the general covariant scalar La- 
grangian L as £ = L{g^i^, g, g^:)y/^. The gauge properties of a field theory under the 
infinitesimal diffeomorphisms ^ = x^ + /S.^x^ , A^x^ = — are expressed through 
a Lie derivative. The latter accounts for the net variation of a field due to both its tensor 
and argument variations. For the metric field it is as follows: 

\9^^u = d^i^gxv + d^^^g^x + i^dxg^.^ = 

= v^e. + v.e^, (1) 

with = g^x^^ and being a covariant derivative with respect to (?^^. It follows hereof 
that 

Aev^ = Iv^g'^'-A^g,, = d^iV^gC). (2) 



2 

By definition, the same should fulfill for /S.^y/—g^. Thus 



A^\Yi^g/g,=i^d^\n^g/g,, (3) 

i.e., g/g^: transforms under G as a scalar. 

Now, restrict the set of S,^ by a subset defined by requirement of invariance of a 
dark modulus (the unimodularity condition)]^ 

A5„v^ = 5M(v^e)=0. (4) 

This singles out in a covariant fashion the residual unimodular invariance subgroup H = 
U. At that 



A^Jn^g = eud,\riJg/g. (5) 



^Restricted ab initio to coordinates where = — 1 (if possible), the covariant unimodularity con- 
dition reduces to the non-covariant transversality condition: = 0, with = ^^|g,=_i. For the 
"transverse/TDiff gravity" based on the latter condition see, e.g., [27]-|29j. 



5 



is not bound to vanish. The general covariance forbids dependence of L on and g 
separately, leaving just 

^ = L{9^u,9/9*)\^- (6) 

An arbitrary g^, is unavoidable in the general covariant framework of the unimodular 
invariant field theory, despite the fact that it might superficially seem irrelevant. Besides, 
it is urgent in practice. In addition to the implicit manifestations it may conceivably result 
in some direct ones. The terms violating general invariance and containing derivatives 
of g/g^, result ultimately in the scalar propagating mode in metric. In the absence of 
derivatives, (7/(7* becomes just an auxiliary field. In what follows the general covariant and 
unimodular invariant metric theory of the bimode/scalar-tensor gravity - the Unimodular 
Bimodal Gravity (UBG) - is consistently exposed]^ 

2.2 Unimodular Bimodal Gravity 

Extended gravity Lagrangian In the framework of the effective field theory any 
choice X = X{g/g^) is a priori allowed as a field variable instead of g/g^,. The differ- 
ent choices are related through the field redefinition. Thus let us take without loss of 
generality as the (dimensionless) scalar-graviton field 

X = lnv^/v/=^. (7) 

This choice proves to be advantageous from the symmetry considerations (see below). In 
this terms, one has for the UBG Lagrangian generically 

-^UBG = Lg + Lg + ALgs + Lm, (8) 

where 

Lg = -^-KlR, L, = ^«:^VX-VX (9) 

present the kinetic terms, respectively, for the tensor and scalar modes, with VX ■ VX = 
g'^'^V ^XV i,X , etc., and V^X_= 9^X. The general invariant Lg is chosen for simplicity 



as in GR (cf. though Sec. 4.2). Here R is the Ricci scalar and Kg = l/y/SnG is the 
mass scale for tensor gravity, with G being the Newton's constant. In the unimodular 
invariant Ls, the parameter Kg is a mass scale for the scalar gravity. For tensor dominance 
of the bimode gravity (generically \Lg\ > \Ls\) one would a priori expect Kg/ Kg < 0{1). 
Moreover, to comply with astronomic observations for galaxies it proves that Kg/ Kg ~ 
10~^ (see Sec. 4.2). Further, ALgg{R, XjV f^X, . . .) is the rest of the gravity Lagrangian 



depending, generally, on R and X. In particular, ALgg may include some scalar potential, 
ALgg = —Vs{X) + . . .. At that, the non- vanishing asymptotically constant part of the 
potential, V^s|oo, may be attributed to the cosmological term. The Lagrangian for gravity 
is to be supplemented by that for the non-gravitational matter, Lm{i', R,X,'V ^X, . . .), 
with ip designating a generic matter field. In principle the latter may correspond both 
to the luminous matter and a putative non-gravitational DM. Finally, instead of being 
imposed explicitly, the restriction ([t]) may be enforced classically with the help of the 
undetermined Lagrange multiplier A by adding to L the term 

^A = A(v^/v^-e-^). (10) 



^To be distinguished from the "plain" (by default, monomode/tensor) UG sometimes used as a 
synonym of UR. Under the phenomenon of unimodular gravity we understand generically both the 
unimodular bimode and monomode/tensor gravities, the latter of them treated as a marginal case. The 
term bimode gravity refines the more sophisticated one - "metagravity" - used previously in [221 123) . 
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This allows one to treat X as an independent variable. 



Global compression symmetry First of all, in the spirit of the effective field theory, 
the extra terms in Lagrangian containing derivatives of X are to be suppressed compared 
to the kinetic term for X by powers of I/k^, and thus may be neglected. To further 
terminate the Lagrangian, enhance the residual unimodular invariance by a dynamical 
global symmetry defined in the fixed coordinates through the field substitutions as follows: 

g{x) g{x) = k^g{kox), 

g*{x) g*{x) = g^{kox), (11) 

with ko > being an arbitrary constant. This is a generic symmetry distinguishing 
the dynamical and non-dynamical fields. For the former ones it coincides (modulo the 
coordinate redefinition) with the conventional global scale symmetry, being a part of the 
general coordinate transformations. For this reason, the general invariant part of Lubg 
{Lg and, supposedly, L^) is global symmetric. A common multiplicative factor appearing 



in £ due to \/—g does not influence the classical equations. Eq. (11) will be referred to 



as the compression transformations. In these terms, the moduli are "incompressible'' 



in contrast to the dynamical metric. Further, it follows from (11) that X transforms 
inhomogeneously under the compressions: 

X{x) ^ X{x) = X{kox) +4:lnko. (12) 

The emergence of the respective (approximate) global symmetry may serve as a natural 
reason for suppression of the derivativeless couplings of the scalar graviton as violating 



the global symmetry (e.g., Vs{X), etc.). Due to (12) the scalar graviton may be treated as 
a (pseudo-)Goldstone boson for a hidden/non-linear realization of the compression sym- 
metry. Associated with the scale invariance, such a scalar graviton/ (pseudo-)Goldstone 
boson may more specifically be termed the s?/sio/onj^ with the tensor graviton being con- 
ventionally just the graviton. The systolon presents a scalar compression mode in metric 
in addition to the transverse-tensor, four- volume preserving deformation mode presented 
by the (massless) graviton. Finally, imposing 

X{x) \{x) = k^\{kox), (13) 

one gets that Lx does not violate global symmetry, too. 

The (approximate) global compression symmetry may be used as the third basic in- 
gredient of UBG, in addition to general covariance and unimodular invariance/relativity. 
To such a reduced version of the theory we will adhere in what follows. This is a minimal 



bimode extension to UR. Namely, from (|9j) and (10) in a formal limit of "switching-off" 
the scalar mode, X — )■ 0, one recovers at the classical level UR, the latter, in turn, being 
classically equivalent to GR with a cosmological term. For more generality, we still re- 
tain the terms Vs{X) and Lm{X) assuming that they are the leading corrections to the 
otherwise global symmetric Lagrangian]^ 

^From the Greek uvaroX'f] meaning the compression, contraction. 

^UBG may be considered as a general covariant counterpart of a theory of gravity [18]-[2Q] comprising 
in metric the (tensor) graviton and dilaton, but restricted ab initio (hereof the term "restricted gravity") 
exclusively to coordinates corresponding in UBG to 5* = —1. The general covariance proves though to 
be crucial for real treating the theory and associating it with DM. In passing, it is the gauge invariance 
which is in fact restricted, the gravity itself being rather extended. 
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Extended gravity equations Varying the total action (assuming ALgg = —Vg) with 
respect to g'^'^ and X independently one gets the tensor and scalar gravity equations, 
respectively, as follows: 



kIV -VX + dVs/dX - Xe-^ = dL^/dX, (14) 

with 

V ■ VX = g'^'V^V.X = d^iy^g'^'d^X)/^ (15) 
being the general covariant d'Alambertian operator. Conventionally, 

T^^, = J^^^'^J^'^^ = IdL^ldg^"" - L^g^, (16) 

is the canonical energy-momentum tensor of the non-gravitational matter, and tg^^ is the 
similar scalar- field tensor given by the unconstrained Lg. 

tsi^v = ^ iJiCy^ - -^o -V^g^y. (17) 

In the above we have conventionally introduced the dimensionfull scalar field a = HgX. 
Introducing Ws as Hg times the canonical inhomogeneous scalar-field wave operator: 

Ws = Kg{V -Va + dVs/da-dLm/da) (18) 

one may present the scalar gravity equation as follows 

Ws = Ae-"/'^% (19) 

with the r.h.s. playing the role of a (virtual) source density. 

Excluding A from (14) we finally arrive at the extended gravity equations in the 
superficially conventional form 



a 



Here T^^y is the total "bare" energy-momentum tensor, with the tensor gravity included 
only in the minimal fashion through metric. At that, the bare effective tensor for systolons 

Ts^y = tsf,y + Wsg^y (21) 

looks like the canonical one, ts^y, except for the Lagrangian potential Vs substituted by 
the effective one 

As = Vs + Wg, (22) 



with Ws given by (18). Otherwise, T^^j, may formally be brought to the form: 

Tsf,y = (ps + Ps)nf,ny + psg^,y, (23) 
where = V ^cr/ (— Vo" • Vo")^/^ [n ■ n = —1) and 

= -V(T • Va/2 + A„ ps = -Wa ■Wa/2- ks- (24) 
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Superficially, this has little to do with a continuous medium, the account for the proper 
gravity contribution being to this end in order (see Sec. 



4.1) 



Finally, varying Lx with respect to A one recovers the constraint ([T]) which is to be 
understood in (20). The non-canonical contribution to Tsfj,u appears ultimately from the 



scalar- field kinetic term due to such a constraint, missing in GR with a scalar field. This 
may be inferred directly without Lagrange multiplier, with Tg^u appearing as a canonical 
energy-momentum tensor from the Lagrangian ^ with the explicit constraint ([T]) 
Presenting (20) in the equivalent form 



R 



■fiu 



,2 \ ^^^•' 



T 



Tq + T/, / = 1, 2, 3, we find in particular that 



(25) 



(26) 



This is of special importance for revealing the energy content of a static space structures 
(see Sec 



reads 



4.1). In the empty space, with = except possibly for a point, Eq. (25) 



R 



(27) 



where is given in general by (22) and (18). Studying this equation and consequences 



thereof in the static spherically symmetric case is the main concern of the following. 
Under the quasi-harmonicity condition, Wg = except possibly for a point, Eq. (27) 



reduces to the conventional Einstein-Klein-Gordon system in GR for a self-interacting 
scalar field minimally coupled to gravity, the solutions for both theories thus coinciding 
in this case. 



Spontaneous global symmetry breaking Due to the contracted Bianchi identity, 
Vfj,{R^ — R/2 6^) = 0, the energy-momentum tensor of the non-gravitational matter and 
systolons collectively, not generally by parts, is bound to covariantly conserve, V ^Tj^ = 0. 
In view of ( 18 ) this results in the third-order differential equation 



d^^Wg + Wgd,a/Kg = -V^T, 



{2i 



which as it stands is hardly of use. However, in the case if Lm is independent of a (in 
particular, in empty space, Lm = 0), the emerging general invariance of the matter action 
implies that V^T^^^ = 0. Thus dy{\n \Ws\ + (J / Hg) = 0, and there appears the first integral 
of motion 

Ws=Aoe-''^^% (29) 



with Aq being an integration constant. It follows then from (19) that now A = Aq. With 
account for (13) this signifies spontaneous breaking of the global symmetry a — ?■ cr + o"o, 
with (To being a constant. In view of (18), Eq. (29) becomes nothing but the bona fide 
second-order scalar-field equation: 



V ■ Va + dAs/da = 0, 



(30) 



with the effective potential looking like the conventional one plus the exponential contri- 
bution 

A, = K + Aoe-'^/'^\ (31) 
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In this case UBG is classically equivalent to GR with a scalar field supplemented by the 
exponential potential in Lagrangian (with an alien for GR scale Ks)j^ At Aq = the 
theories moreover classically coincide. For physical reason, there should fulfil Aq < (see 



Sec. 4.1). This means that due to a potential well produced by a systolon condensate, 
the local vacuum with the spontaneously broken global symmetry lies lower than the 
symmetric one (Aq = 0), ensuring ultimately spontaneous breaking of the symmetry. 

UBG and beyond Let us make some remarks posing UBG among the related metric 
theories of gravity with the additional degree (s) of freedom, which may have bearing 
to DM. 

(i) Unimodular scalar-tensor gravity First of all, putting g^^ = {g / g^Y^^gunv, with 
det^fu^,, = gu = g*, one may equivalently present (|6]) as £ = L{gu^v, X)y/—g^, where 
X = X(g/g^) may be chosen, e.g., as in ([T]). Under the global compression symmetry one 
has gu^ui.^) Qu^iuix) = gufiuikox), while g^: and X transform as before. Now choosing 
a = KgX as an independent scalar field one may reduce the original theory to the 10-field 
unimodular scalar-tensor gravity, which proves to be nothing but the unimodular bimode 
gravity in disguise. The inverse transition is achieved by putting gu^^, = e''^^'^'^" g^u and, 
respectively, a / Hs = (1/2) In (7/ (7=,,. The two theories may thus be treated just as the 
generic UBG in the two different frames: the bimode and scalar-tensor onesj^ 

(ii) General scalar-tensor gravity More generally, one may confront UBG, contain- 
ing 10 dynamical fields, with the general invariant theory defined by the arbitrary C = 
L{gfiu, o')\/—g containing additionally an independent scalar field a, altogether 11 dy- 
namical fields, ^^P^ As it was stated earlier, under the separate conservation of the ener- 



gy-momentum tensor of the non-gravitational matter (in particular, in the matterless 
vacuum) UBG is classically equivalent to the general scalar-tensor gravity, with a sup- 
plementary exponential contribution to the Lagrangian potential of the latter emulating 
the proper spontaneously emerging term in the classical equations of the former. Never- 
theless, in the context of DM there are two important differences. First, to match with 



astronomic observations the sign of the extra contribution is to be negative (see Sec. 4.2). 
Being harmless in the classical equations such an unbounded from below term in the La- 
grangian could result in a quantum inconsistency. Second, appearing as an integration 
constant in the classical equations, Aq may vary for different solutions, whereas in the 
Lagrangian potential it should be fixed ad hoc ones forever. Thus in the DM context 
UBG seems to be more safe and fiexible compared to the general scalar-tensor gravity. 

(Hi) Bimetric multimode gravity In a wider perspective of the general relativity vi- 
olation one may consider the multimode gravity, with the general covariant Lagrangian 
density £ = L{g^^, g^fj_^)y/^ depending on the dynamical and non-dynamical metrics 
g^u and g*^u, respectively, as well as their determinants. Basically, such a theory is char- 
acterized by its residual gauge invariance H C G. Containing 10 independent dynamical 
fields, the theory with the trivial H = I may thus encounter up to 10 physical degrees of 
freedom (including, in general, the ghost ones) and present potentially a lot of problems 

®This statement resembles a similar one about the classical equivalence between UR and GR supple- 
mented by the cosmological term in Lagrangian. This can be seen at Vg = in a formal limit cr — > 0. 
^In the latter frame the unimodularity literally means the invariance of gu- 
^°For a confrontation between the transverse/TDiff gravity and GR with a scalar field see [30] , 
^^As a particular implementation of such a theory there may be mentioned the well-known Brans- 
Dicke model [31]. In the Einstein frame the scalar field in this version of the theory acquires the direct 
derivativeless coupling with matter being thus strongly restricted observationally (see, e.g., 
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and ambiguities. To eliminate/reduce them the residual unimodular invariance, H = U, 
may, e.g., be imposed. In conjunction with general covariance this would retain the de- 
pendence of L just on g^i^ and g/g^,, with g^: = det((7*^j,), and no more than a single extra 
graviton. Nevertheless, in the wake of the developed approach to DM one could envisage 
in this direction a multimode gravitational DM,^^ 



3 Equations and solutions 

3.1 Static spherical symmetry 

Quasi-Galilean coordinates Consider the static spherically symmetric field configu- 
rations. Start with the quasi- Galilean coordinates = {x^,x"^) = (t,x) (m = 1,2,3), 
where the general covariant line element is as follows: 

rfs^ = adt^-(6-c)(nrfx)2-crfxl (32) 

Here n (n^ = 1) is given by n"^ = x"^/\x.\, with |xp = x^ = 6mnX"^x"'. The metric 
variables a, b and c are arbitrary functions of |x| alone. The same concerns g^, (and 
thus a). The respective metric looks like 

900 = fl'mn = ~bnm^n c{6mn ~ ^m^n)- (33) 

where for short we designated Um = Smn^"' (vs. conventional = gmnn"')- The rest of 
the metric elements is zero. Respectively, the inverse metric is 

= ^, ^'"'^ = nV - ^(5"^" - nV). (34) 

Due to the rotation invariance one can choose the spatial coordinates in a point x 
so that n = (1,0,0), bringing the metric in this point to the diagonal form (yf^j^ = 
diag (a, —6, — c, — c). Thus ^/—g = Vabc, and the scalar field looks like 

a = h^s\n{^/abc/^y—g^:). (35) 

The quasi-Galilean coordinates are appropriate for revealing the energy content of the 
static space structures (see Sec. 4.1). 



Polar coordinates For real dealing with the spherical symmetry better suit the polar 
coordinates x'^ = {x^,x"^) = (t,r,9,ip), m = r,9,ip (with a unit of length /q tacitly un- 
derstood where it is necessary). Conventionally, x = rn = r{sm6 cos(f, sm6 simp, cos6). 
The line element now reads 

ds"^ = adt^ - bdr^ - cr^dQ^, dn^ = dO'^ + sin^ ^c^(^^ (36) 

with the metric being g^i^ = diag(a, —b, —cr^, —cr^ sin^ 9). The scalar field looks as before, 
with the transformation Jacobian cancelled out in the ratio of the two scalar densities of 
the same weight. The functions a, b, c and g^, depend on r alone. 

""^^For the bimetric gravity in the context of the massive tensor graviton, with a priori unspecified 5*^1/ 
and extra terms only in the Lagrangian potential, cf., e.g., j33l I34j . 
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The systolon bare effective energy-momentum tensor (17), (21 ) and (22) is then Tg'^ = 
diag(po, -Pr, -Pe, -p^), where 



Po 



-Pv 

Pr 



26 
1 

26 



(37) 



r, 6', the bare pressure and prime 
indicating a radial derivative, There follows hereof, in particular, an important spatial 



with Po being the bare energy density, p„ 



invariant (see Sec. 4.1): 



rp \) rp n 

''-SO J-Sn 



Po + Sp„ = -2A,. (38) 

With being diagonal and = Rg, there is left three independent gravity equations 
which may be brought to the form: 



Rq 



Rq 



Kq - Kg 



-^(t ° — T " — 2A 

— (t 0- 

— (t °- 



J-mr + 



(39) 



In the above the curvature elements are expressed through a, b and c. The contracted 
Bianchi identity reads 

W: + W^a'/K, = -T'j;. (40) 

Under the general invariance (cx-independence) of Lm (in particular, at Lm = 0) one 
has T'j^l = 0, and there appears the first integral of motion Ws = Aoe~'^/'*% with Aq 
being an integration constant. The third-order equation (40) reduces in this case to the 
second-order scalar-field equation (30) which looks now like 

V ■ Vex = —( \/a/bcr^a'^ / Vabcr^ = 

(41) 



^a/bcr'^aj 
(Ao/«.)e-'^/''» - dVs/da. 



Radial coordinate/gauge fixing The choice of the radial coordinate r is not unique. 
Under the local radial rescaling r — )■ f = f(r), the general covariant line element is 
moreover formal invariant: ds"^ = ds^ = ddt^ — bdr'^ — cf'^dQ'^, with the relation 

a{r) = a(f(r)), 

6(r) = {dr /dryb{r{r)), 

c(r) = {r/rYc{f{r)). (42) 
The quasi- Galilean ^/—g = \/abc transforms under rescaling as 

V^=\/=^^f (43) 

(and so does y/—g*), with cr transforming as a scalar, a{r) = a{f{r)). The three- volume 
element y/—gd^x = y/abcr'^drdil'^, being (spatial) invariant, is also formal invariant. It 

""^^Here and in what follows, we define for simplicity the energy densities without y/~g. 
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follows from (42) that metric, asymptotic Minkowskian with respect to r, retains this 
property with respect to f if there fulfills asymptotically f = r(l + 0{l/r)). Thus, the 



ansatz (36) does not allow to fix the metric uniquely. 



The aforesaid ambiguity can be removed by the radial coordinate/gauge fixing. Namely, 
the three gravity equations (39) contain superficially four variables: a, b, c and a, only 



three of them being independent due to constraint (35). To account for the latter, two 



opposite routes of dealing with a priori unknown ^f* are envisaged. 

{i) Explicit The direct route is to assume in some starting coordinates and solve 
the equations with the proper boundary conditions. Then postulating a relation between 
the starting and observer's coordinates, transform the solution to the latter coordinates, 
including boundary conditions. In practice, one may start from the transverse coordinates 
defined hy g^ = —1, what is equivalent to imposing the scalar-field dependent gauge 
y/—g = . In the arbitrary observer's coordinates the solution will thus explicitly 
comprise g^: ^ —1. 

(m) Implicit g^ The inverse route is to get (virtually) rid of g^ in favour of a. Imposing 
then a gauge F(a, 6, c) = as a definition of observer's coordinates, solve the equations, 
with the proper boundary conditions imposed already in the latter coordinates. The 
dependence on g^ remains implicit. Having found the solution one can extract hereof 
the required g^ solving, in a sense, an inverse problem. Inserted back into the gravity 
equations this is bound to produce in the direct route in the same observer's coordinates 
precisely the given solution. Having a priori no knowledge about g^. we will adopt in what 



follows such an inverse route, with several particular gauges compared (see Sec. 3.2) 



Reciprocal gauge It proves to be convenient to impose the gauge ah = \ (refer to 
it as the reciprocal one) Note that as it stands this gauge is not invariant under the 
dynamical global symmetry (11). Designate A = a = 1/6, C = r^c (some element of 
length Iq is put to unity here). 



With account for 



1 [CA)' 

2 C '■ 

1 {CA)' 

2 C 

1 {AC')' 

2 C 



A 

1 
C 



C" 

~c 



\9l 

2~C^ 



(44) 



and V ■ VX = —{ACX')'/C, the gravity equations look like 



CA' + CA 



(InC)' 



CA" - AC" 



v'Mcx'Y - 



+ 



c 



rp rp n 

-L niQ -L nin } i 



2 kIA^"'^ 



-2 



2C/ o_ 



(45) 
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The ultimate reason is the appearance in this gauge of an exphcit exact solution valid also in GR 



with a scalar field (see Sec. 3.2 ) 
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where v"^ = 2k1/k'^. With Hg taken as an overall mass scale, Vg is the single free parameter 
of the vacuum gravity Lagrangian (under Vs = 0). According to (45) v"^ plays the role 



of a coupling between the scalar and tensor gravities. (For its observational meaning see 
Sec. 4.2 ) With 1;^ ~ 10"'^ the coupling proves to be weak. 



Lagrange multiplier ansatze In the equations above, the Lagrange multiplier is ex- 



cluded. It may be revealed through (19) as follows: 

- uliACXj/C + dVs/dX - dL^/dX = Ae"^. (46) 



Given a solution to (45 ), Eq. (|46|) defines A. V.v., making an ansatz for A (to be confirmed) 



one can look for a respective solution to (45) (if any). This allows one to (partly) 



disentangle the scalar and tensor gravity equations. 

Being interested in the case of the empty, but possibly for the origin, space we envisage 
at Lm = the three following ansatze. 

(i) Singular ansatz Here A = Aq, with Aq being a 5- type function concentrated at 
the origin and determined implicitly through self-consistency. 

(ii) Regular ansatz Here A = Aq everywhere (including the origin), with Aq being an 
arbitrary constant. 

(iii) Interpolating ansatz Here 

A^/fo, r<ro, 

I Ao, r > ro, ^ ^ 

with ro being some matching distance. In reality A should be smoothed around tq. 
Consider these ansatze in turn (neglecting by V^). 

3.2 Quasi-harmonic solution 

Reciprocal gauge The singular ansatz A = Aq reduces to the quasi-harmonicity con- 
dition V ■ VX = 0, or [ACX'Y = (except for the origin). The exact solution to the 
proper equations was found previously in [24]. Here we reproduce it in the nutshell. 



Combining the first and the last gravity equations (45) one gets (AC)" = 2 with the 
ensuing relation 

AC = A = {r -ri){r -r2), (48) 

where the roots vi and r2 are a priori either real or complex-conjugate of each other. 
Introduce a new radial coordinate 

X^-/f, (49) 

with X = 1/r + 0{l/r'^) at r — )• 00 (a unit of length /q is understood here and in what 



follows). In these terms the first gravity equation (45) and the scalar-field equation (46) 
reduce to the free harmonic form 

In A/dx^ = d^X/dx^ = 0, (50) 

with general solution 

In A = In Ao - uox, X = Xq - qoX, (51) 

depending on the integration constants Aq, Xq, z/q and ^q. Conventionally, impose the 
asymptotic free boundary conditions A = 1 and X = Oatx = 0(r— 00), so that Aq = 1 
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and Xo = 0. Substituting InC = In A — In A into the second gravity equation (45) we 
finally find the integration condition 



It follows that for solution to be real it is necessary that roots be real ((ri — r2) > 0) 
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In this case we get from (49) at ri 7^ r2 (let by default ri > r2): 



1 , r — ri , , 

X = In 53 

Ti — r2 r — r2 

(modulo a constant), or inversely 

ri + r2 ri-r2 ^, (ri - r2)x 

r = cth . 54 

2 2 2 ^ ^ 

At ri = r2 = vq one has x = ~ ''^0) and, respectively, r = vq + 1/x- Altogether, the 
respective line element and the scalar field look like 

X = ^Ing, (55) 

where 

g = (r - ri)/(r - r2), (56) 
with u = 1^0 /{r I - r2), = ?"o/('"i - ^^2) and z/^ + t;^^^ = 1. 



The gravity and scalar-field equations in the reciprocal gauge, (45) and (46), still 
possess in the vacuum a residual invariance under the global shifts of r. Due to the 
invariance one can redefine r as r — )■ r + r2, so that q = 1 — r f /r and A = r^g, with 
r/ = ri — r2 > 0. The looked-for solution (designate it with a subscript /) acquires the 
following standard form: 

a; = 1/6; = {l-Vf/ry^, 

Cf = (l-rf/ry-'^f, 



V2af/K, = J:j=VsXf = rij^l-u]\n{l-rj/r), (57) 

with Uf = Uo/vf and a signature factor rjf = ±1. The solution depends on two canonical 
parameters rj and i^/.^^ 

For consistency there should fulfil |z//| < 1. The solution is unique (up to constants 
Aq and Xo) and bound to be singular, with the regular solution being necessarily trivial 

{vf = 0). Note that the particular case rjj = 1, uj = Vs/^Jl + results in the relation 
In a/ = VgXf, which proves to correspond to a matterless static space structure (see 



Sec. 4.1). The two-parameter solution above supersedes the one-parameter Schwarzschild 
solution for BHs (z// = 1, Xf = 0). In GR it describes BHs dressed with a (free massless) 
scalar fieldj^ Reflecting a singularity in space, with the coherent scalar field treated as 

-'^^In the case of complex-conjugate roots ((ri — r2)^ < 0), the reahty of metric (fg > 0) imphes <;o < 
meaning a ghost scalar-field solution {X'^ < 0). 

-'^^A more symmetric form of the solution (still in the gauge a5 = 1) would correspond to the shift r — > 
r+{ri+r2)/2, so that x = -(1/r/) ln(r-r//2)/(r + r//2), r = (r//2)cth (r/x/2) and A = r'^-{rf/2f. 

^''in GR the respective solution was first obtained in an implicit form in the gauge c = 1 by Bergmann 



and Leipnik [3S] (see also [Mj)- In the explicit form ( [57| it was discovered somewhat later in a different 
context by Buchdahl [37] and was then extensively studied in literature including modification by the 
scalar-field self- interaction (see, e.g., [38]-|45)). 
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DM, such a scalar-modified BH may conveniently be referred to as the dark fracture. 
Accounting for (35) one finds the respective dark modulus as follows: 



-g^f = {l-rf/ry-^f, (58) 
where 

^f = Uf + T]f^l-uj/Vs. (59) 

Now, it is possible to reverse the line of reasoning by saying that the dark fracture is the 



static space structure which originates from the dark modulus (58) and is given by (57), 
with Uf = Uf{nf). In particular, to the constant modulus (/i/ = 1) there corresponds BH 
(z// = 1). At -C 1 there fulfills i^/ ~ 1 — v1{^f — 1)^/2, with Vf being close to unity in 
a wide range of /i/. Thus, the appearance of quasi-BHs (|z// — 1| ^ 1) is parametrically 
enhanced. For an arbitrary fracture the modulus, being complex singular, can not be 
brought to g^f = — 1 by any real coordinate transformations. 

In general, the solution is complex interior to ri and, as such, should be treated 
here as an analytical continuation from the exterior region r > ri into the complex r- 
plane with the cut (— cxd, ri). Moreover, the event horizon at r = ri proves to be point-like 
singular [38]. To avoid the complexity one could perform the admitted shift r — )■ r + ri, so 
that q = 1/(1 + rj/r) and A = r^/g. At r > the solution becomes now real and regular, 
with a singularity only at r = 0. The same concerns the modulus, which may now be 
brought to unity (but for r = 0). The asymptotic of the solution being unchanged under 
the shift, the latter may be considered as a kind of "realization". Being mathematically 
equivalent, such a "truncated" real form seems physically more reasonable. Nevertheless, 
ultimate treating the complex singularity would be of considerable interest 



18 



Until stated otherwise, we pursue the standard form of the solution. Decomposing 
the latter in 1/r one gets in the leading approximation: 

Cf = l-Tc/r, 
T.f=VsXf = -rjr, (60) 



where Vc = ^Jr"^ + ri — Vg. Here and are two independent phenomenological param- 
eters 

^^9 = ^' '^s = Tlf^^-^}rf, (61) 

or inversely 

Taking as independent rf and one has 

p^ = rg/rf, rc = rf-rg, r^ = Vf^r'j - r^. (63) 

For definiteness, we consider the case > > 0, with remaining sign-indefinite 
[rjf = ±1). The parameters and fix the Newtonian-Coulombic approximation and 
have the meaning, respectively, of the gravitational and scalar radii of a fracture. Under 
the radial rescaling r — )■ f(r), with f = r(l -|- 0{l/r)) asymptotically, these parameters 



-'^^For treating the interior of BH as unattainable due to its event horizon shrinking to a point singularity 
by dressing with a massless scalar field minimally coupled to gravity cf. [55) . 
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are invariant and may serve as a substitute to the canonical ones, rj and Vf, with rj being 
the radius of a fracture and Uf its "nudity". (BHs, Uf = 1, are the "nude" fractures.) 
Now we may reverse the approach. Namely, in an exterior region where = let us 



look for an asymptotic free solution to (45), with Vg = and V- VX = —{ACX')'/C = 0. 



With E = VgX, decompose the solution in a power series of 1/r. Starting from (60) with 
a priori arbitrary parameters Vg, Tc and r^, we get in the second order the restriction 
rc{rc+2rg) = r^. Afterwards we can step-by-step uniquely reconstruct the quasi-harmonic 
solution (57) in the empty space within the series convergence region r > rj. Because 



we have nowhere used any assumptions about the distribution of matter interior to r^, 
but for spherical symmetry, the exterior vacuum solution is to be uniquely determined 
by the two interior integral characteristics corresponding to Vg and r^, independent of 
the details of the interior distribution. Physically, Vg reflects the net gravitating energy 



of a fracture and its net systolon energy (see Sec. 4.1). If some "reasonable" matter 
distributions exist, but for the point-like one, this would extend to UBG the Birkhoff 
uniqueness theorem in GR. The Buchdahl solution may thus be to dark fractures as the 
Schwarzschild solution is to BHs. 

In the above we have in fact shown that the gauge most appropriate to the problem 
at hand from the point of view of its physical content is the harmonic gauge. It may be 



given by requirement V ■ VX ~ d^X/df^, or in view of kn by {a/by/^cf^ = 1, with f 
substituting x- (Pot completeness though, we should have found b and c.) The original 
reciprocal gauge (fe = 1/a) may, in turn, be more appropriate for studying the quasi- 
harmonic solution in the exterior region. Below we shortly compare the exterior solution 
in several radial coordinates/gauges peculiar from supplementary points of view. 

Astronomic gauge The coordinates where c = 1 are peculiar in astronomy, with the 
surface of a concentric sphere being conventionally S = Airf"^. In view of (42) such a 
gauge results in the exterior coordinate transformation 



f = \J c{r) r = r(l — rj/rY^ "f)/"^^ r > rf. (64) 

The inverse to this relation being given only implicitly, the exact solution in these coor- 
dinates can not, unfortunately, be presented in an explicit form (but for z/j = 1). 

Unimodular gauge The gauge \/—g = {abY^'^c = 1 is peculiar by the fact that here 
the modulus is directly reflected by scalar field, y/—g* = e~^. In view of (43) this gauge 
results in the exterior coordinate transformation 

r = (s J c{r)r'^dr^ ^ = 
= (3 y"(l - Tf/rf-^fr'^dr^ r > r/, (65) 

with the integration constant set to zero to ensure f = r^l + 0{rf/r)^ at r > Vf. 

Transverse gauge This gauge is peculiar from the theoretical considerations. It is 
given by ^=1, = —1, with the group of diffeomorphisms being the transverse one {d^C,'^ = 0) 
in the respective coordinates. Otherwise, this is equivalent to the gauge y/—g = 
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(which may also be imposed in GR with a scalar field). Applying (43) to g^, one gets the 
equation 

= \/=^=V^^"''^^^ (66) 

so that in the exterior region 

= (3 y"(l - rf/ry-f^fr^dr^ '\ r > rj. (67) 

with the asymptotic relation f = r{l + 0{r f / r)^ . In general, the exact solution can not 
be presented in this case in the explicit form (including GR with a scalar field) . For BHs 
{pf = Hf = 1) all the aforementioned gauges/coordinates (but for harmonic) coincide 
identically implying in particular = —1. 

Isotropic gauge This gauge is peculiar from the phenomenological considerations. In 



view of (|33|), under the gauge c = b one has in the quasi-Galilean coordinates gmn = —b6„ 
fm, n 



1, 2, 3), with ds"^ = adt^ — bdic^. Such a conformal fiat spatial metric is natural in 



confronting with Newton's dynamics. With account for (42) this results in the exterior 
coordinate transformation 

r /b{r)\^/'^dr r 1 dr , , 



so that 

or inversely 

This gives 



(1 + ^Jl-Tf/r^ , r>rf, (69) 
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r = r(l + ^y, f>rf/A. (70) 



t, = v,X, = 2,„v^ln((l-^)/(l + ^)). (71) 

with being the odd function oir/vf. Now the solution preserves the dynamical global 
symmetry, as well as possesses additionally a hidden symmetry: — >■ — tj, z/j — )■ —Vf 
and Tjf — )■ —rif. The same concerns g^. 

Post-Newtonian approximation Without loss of generality present the fracture so- 
lution in the isotropic polar coordinates (with the hat-sign omitted) as follows: 



a/ 


= 1-^ + 

r 


^/ = c/ 


T 

= 1 + 7/- 
r 







2 

9 

2 r2' 



(72) 
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Here /3j and 7/ are the effective parameters for a dark fracture in tlie parametrized post- 
Newtonian (PPN) formalism. It may be said in general tliat(/3 — 1) reflects empirically 
the degree of non-linearity in the superposition law for gravity, while (7 — 1) corresponds 
to the amount of space curvature produced by the unit rest mass |32]. The parameter ^/ 
may be treated as an effective form-factor in the scalar-field Coulomb law. Decomposing 



(71) in 1/r and using (62) one gets in the leading approximation 



with 



/3/ 
If 



1 + 



(1^ 




3rg 




9 ^ 


8 r 




A) 


?>rg 


3 ^ 


8 r 


1 + 






4J 
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Ts/rg = r]f^V^\ 



1. 



(73) 



(74) 



At small e/ a dark fracture closely reproduces BH (r^ = 0, i// = 1, e/ = 0), with the 
scalar dressing being thus very neatly hidden. The precision local tests of GR in Solar 
System result typically in |/3 — 1| < 10"^ and I7 — 1| < 10~^ [32]. With a value near the 
Earth rs/r c:^ 2 x 10~^, rs being the Sun gravitational radius, this implies for the Sun 
as a dark fracture just \es\ < 10^, or i>s > 10~^, the Sun scalar form-factor ^5 remaining 
practically unity. At face value this gives very loose restriction,™ Other observational 



manifestations are conceivably needed to test the theory (see Sec. 4.2) 
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3.3 Regular non-harmonic solution 



Vacuum scaled distance To begin with, take InA = a and Inc = ( as the new 



metric variables and present the gravity and the scalar-field equations, respectively, (45) 



and (46) in the empty, but possibly for the origin, space equivalently as follows: 



^/ a+C^2^/ _ 2^ 

dry dr^ 



dr\ drj 2\dr) 
d f^+^^2 d{a-0 
dr \ dr 

dr \ dr J 



_vl( d_xy 

2\ dr) ' 



2r(l 

( 



K 



OK 



^AdX 



Xe 



-X 



(75) 



It becomes now evident that the second-order system above always has one first integrals, 
with one more appearing at Vg = 0. Moreover, in the latter case the first equation above 
may always be satisfied with a = v^X (modulo a constant), with the set of solutions to 



the remaining equations being definitely not empty (cf., e.g.. Sec. 3.2). Such a solution 
proves to correspond to a matterless static space structure (see Sec. 4.1). 



-'^^The validity of decomposition in 1/r implies that at the distances at hand there should take place 
e/l < (r/rg)i/2 ~ 10^ being still safely fulfilled. 
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For discussion of the transverse gravity vs. observations cf. also |28j . 



19 



Now, under the regular ansatz A = Aq < everywhere (including the origin) let us 
introduce in the empty space a characteristic length scale Rq through 



- Ao = Gk^JRI = ?,{vsKjRof (76) 
and choose the scaled distance r = r / Rq as an independent variable. In these terms the 



vacuum scalar- field equation ( 75 ) at \4 = becomes 



d 



drx dr J 



(77) 



while the tensor gravity ones remain unchanged modulo substitution r — )• r. 



Approximate equations Of particular phenomenological interest is the case t;^ <^ 1 
(see Sec. 4.2). Assuming in this case X = 0{1) and |a|, \(\ <^ 1 (to be confirmed) present 



the scalar-field and tensor gravity equations in the linear in a and ( approximation, 
respectively, as follows: 



d / 2d'X\ 
drK dr J 

d^f2(K 
dr \ dr 



with the last tensor gravity equation 

d / 2 d , 
dr \ dr 



6r^e 



V 



2 / dxy 
~d^j ' 



- 2r(a + C) 







(78) 



(79) 



serving as a consistency condition. Clearly, the coupling of the tensor and scalar gravity 
modes is weak, whereas the self-coupling of the scalar mode proves to be strong. The 
second equation (78) has the general solution a — v^X = ci/r + cq, with ci and Cq being 
some constants. 

The driving equation in the system above is that for X. Having solved the equation. 



one can easily find a and C- It follows from (78) that if X(t) is a particular solution 



then an equivalence class of solutions may be obtained by the inhomogeneous scaling 
transformations: 

X{t) X{t) = X{koT) - 2 In fco, (80) 

with fco > being a constant. This is a reminiscence of the dynamical global sym- 
metry (12) in the global non-symmetric gauge ab = 1. At face value, (80) reduces to 



reparametrization Rq — Ro/ko supplemented by a shift in X. Because i?o is arbitrary 
and a shift in X does not matter at \4 = 0, we restrict ourselves by fco = 1- 
To study the scalar-field equation, introduce the new variables 



t = In 3t^ 



X-t 



and present the equation equivalently as follows: 



d'^Z 1 dZ 



1 

2'^ 



^1) 



^2) 
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Figure 1: Phase plane {Z, Z): arrows designate the normahzed direction field determined 
by dZ/dZ. The trajectories {Z{t),Z{t)) are tangential to the direction field everywhere. 
Solid line - the regular trajectory (Z — )■ — 1 at t — )■ — cxd); dashed lines - irregular 
trajectories {Z — )■ — oo at t — )• — oo). The origin Z = Z = is attractor at t — )• oo. 

Putting further Z = dZ/dt one can bring the second-order ordinary differential equation 
above to the autonomous first-order differential system: 

^ = Z 
dt 

dZ 1 • 1 / 7 \ , , 

Such systems are known to be basically characterized by the types of their exceptional 
points given by the requirement dZ/dt = dZ/dt = 0. The system above has the single 
such point, Z = Z = 0, and the latter proves to belong to the stable focus type. The 
respective phase plane (Z, Z) is presented in Figure [I] It clearly shows the attraction 
point and a distinguished trajectory (the solid line) to which all other trajectories (the 
dashed lines) tend to accumulate departing it nevertheless sooner or later. 



Accordingly, there are three classes of solutions to ( 78 ) : 

(z) an exceptional solution reflected by the attraction point at the origin of phase 
plane; 

[ii) a regular at the origin (r = 0) solution corresponding to the regular trajectory 
(Z — > — 1 at t — )■ — oo); 

{Hi) the irregular at r = solutions corresponding to the irregular trajectories {Z — ?■ 
— oo at t — )■ — oo). Consider them in turn. 
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Exceptional solution The exact exceptional solution corresponding to Z = Z = 
looks like 



X 



IK 



In 3r2, 
-InSr^ 



with the integration constants properly chosen. This solution may serve as a reference 
one, with all other solutions approaching it at r — )■ cxd. Note that X is invariant under 



the rescaling (80). Inverting (35) one gets in the leading i;s-order the hard-core dark 



modulus corresponding to X: 



-X 



l/(3r2). 



Implementing (43) to g^, gives for the transverse coordinates the cuspy relation f 



(85) 

-1/3 



Regular solution Decomposing a regular solution in the powers of r (only even powers 
prove to enter) we get the looked-for solution (endow it with the subscript K) as follows: 



and 



Xh = ah/vl 



Ch/vI 




10 ^30 



r< 1, 
r > 1 



10 ^35 

In 3r2 + 2, 



r< 1, 
r > 1, 



^6) 



^7) 



with X/i(0) = being imposed. The restriction (79) fulfils identically with the same 



accuracy both at r < 1 and t ^ 1. Note that though Xh changes under the global trans- 
formations ( [80| its asymptotic remains invariant. The solution satisfies the matterless 
condition, ah = v^Xh (see Sec. 4.1). The respective static space structure is nothing but 
the dark /ia/oj^ The parameter Rq plays the role of the soft-core radius (see Sec. 4.2). 
The equations above can be extended analytically to < (corresponding to Rq < 0, 
Aq > 0), though with loosing the halo-type solution. This explains the earlier made 
choice Aq < 0. 



Inverting (35) one gets in the wide region of r, where = ^{'^s) 1' soft-core 



dark modulus as follows: 



-X, ^ fl-r2 + 0(r4), r<l, 
\l/(3r2), r»l. 



Implementing (43) to g^,h gives then for the transverse coordinates the relation 



-1/3 



r<l, 

r > 1. 



(88) 



^9) 



To study the asymptotic behaviour of X^ in more detail put X^ = X + AXh- In these 
terms the approximate scalar-field equation looks equivalently like 



d / 2d^Xh 
dr \ dr 



(90) 



^^Indeed, equating the Newtonian gravitational attraction force Fg — ma'f^/2, acting in the given 
metric on a test particle with the mass m, to the centripetal force Fc — mvf^/r, corresponding to the 
circular rotation velocity Vh-, we could already anticipate the asymptotic constant = Vg, characteristic 



of the galaxy dark halos (see Sec. 4.2). 



22 




Figure 2: Regular solution Xh. bold solid line - numerical result; dashed line - piece- 
wise analytical approximation with 60 and tq as in Figure 3. Thin line - exceptional 
solution. 



Assuming |AXft| < 1 (to be confirmed) and retaining the linear in AXh part we get the 
solution at r > 1 as follows: 



AXh = (So/ a/t) cos {{V7/2) Inr/ro 



(91) 



Here Sq and Tq are some integration constants to be inferred from matching with solution 
at r < 1 or from comparison with the numerical solution. Clearly, X^ oscillates with 
attenuation around X approaching the latter at r — )■ 00. The regular solution is unique 
and may be prolonged in the arbitrary order in to the solution of the exact equa- 
tions (75). The behaviour of Xh is shown in Figures |2| and [31 It is seen, in particular. 



that Xh gets strong already at the moderate r, so that the account for potential Vs{X) 
may become important in this region,^ 



3.4 Irregular non- harmonic solution 



At last, consider the interpolating ansatz for A given by (47). To satisfy this ansatz. 



drop off the requirement of regularity at the origin and look for a solution interpolating 
between the two previous solutions. Put without loss of generality X = Xh + AXh- In 
these terms the approximate scalar- field equation corresponding to A = Aq at r > ro/-Ro 
looks equivalently like 



d / 2d AXh 
dr \ dr 



6r e~ 



Similarly, putting a = ah + Aah and ( 
equations at r > ro/-Ro as follows: 



^'^(e-^^'^ - Ij. (92) 
h A(h present the approximate gravity 



d 



d 



dr V dr 



(Aah-vlAXh)^ = 0, 



d f^dACh 



dr V dr 
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For additional details on the regular solution see [25l [26 



23 



-0.6 



Figure 3: The asymptotic of the regular solution Xh- solid line - numerical result for 
AXh = Xh — X] dashed line - analytical approximation with 5q = —0.75 and tq = 2.0 

2 ^fdXh , l dAXh \ dAXh 
~ ""'^ [ dr ^ 2 dr ) dr ' ^^^^ 

with the consistency condition remaining as before: 



d / 2 d 
dr \ dr 



{Aan - ACh) - 2r(A«ft + A^)) = 0. (94) 



Assuming |AX;i| < 1 and \dAXh/dT\ < \dXh/dT\ (to be confirmed) consider the linear 
in AXh approximation. Decomposing AXh in the series in r starting (by assumption to 



be verified) from 1/r and accounting for (86) we get the solution as follows: 



~ \ (So/V^) COS {{V7/2) In r/ro) , r > 1. ^^^> 

In the above, Si is a small normalization parameter to be fixed by subsequent matching 
with the dark fracture solution. Likewise, 6o and tq are some integration constants which 
can, in principle, be fixed by further applying the perturbation procedure in AXh- These 



constants can be estimated by matching the two branches of (95) or from comparison 



with the numerical results. Solving (93) we then get 



Aah = -- + v^,AXh, (96) 

T 

At - _k_„2jr-lr' + 0{T'), ro/Ro<T<l, 
" r ^^'^\0(l/v/^), r>l, 

with 6a and 6c being some small parameters to be fixed by further approximation. The 



consistency condition (94) is fulfilled with the same accuracy. Several representative 
solutions corresponding to the irregular trajectories somewhat close to the regular one 
are shown in Figure El The less \6i\, the better is the approximation, 
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^^The decaying behaviour of the irregular solution with (5^ > is superseded eventually by the growing 
one at the tiny r (not shown), what lies though beyond the region of approximation. 
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ffj = -0.12 
5j = -0.07 



IX 0.8 1.0 



<5j=0.0fi 

i5i = 0.12 



Figure 4: Representative irregular solutions X: solid lines - numerical results for AXh = 
X — Xh] dashed lines - analytical approximation with the respective 6i for the pole 
singularity. 



To arrive at a solution valid approximately in the whole interval of r, match the 1/r- 
pole term in AXh with the respective term in Xf, the latter corresponding to A = Aq at 
r < tq/Rq (and similarly for Aah and Aa^). In view of (60) this implies 



5i 


= rj{v 






= (rg- 






= (Vi 


' s 



Rn 



(97) 



In particular, the case Sa = {j-g = VsTg) proves to correspond to a matterless fracture, 
with In a/ = v^Xf (see Sec. 4.1). Partite further AXh as AXh = Xf + Xt, where 
Xf = —6i/T = —rs/{vsr) is the scalar tail of fracture and the rest, Xj, is attributed to 
fracture-halo interference. Altogether, the total solution (designate it by the subscript /) 
can be presented as a coherent sum of three contributions, Xi = Xh + AXh = + + 
Xi = Xf + Xheff (and similarly for the metric components ai and (i). At r — oo the 
interference disappears due to disappearance of Xf. At 5j = (r^ = or i?o oo) it 
disappears identically. The behaviour of Xi/5i is show in Figure |5} Finally note that in 
the Laurent decomposition of the exact X/, the part comprising powers of 1/r may be 
associated with Xj, the even powers of r with Xh, and the odd powers of r with Xj. 

The gravitational potential ai reveals the property of confinement, with the reciprocal 
behaviour, —Tg/r, at the tail of fracture superseded eventually at the periphery of halo by 
the logarithmic potential of gravitational attraction, t;^ InSr^/i^^. The confinement scale 
Rq is otherwise the halo soft-core radius (see Sec 4.1). Call the respective static space 
structure the dark lacuna,"^ Due to the coherent scalar field getting strong towards the 
periphery of halo, the gravitational confinement within a lacuna should in reality be only 
partial, being terminated ultimately by the potential Vg and/or the influence of the near- 
by lacunas. Generally, the lacunas depend on three distance scales which may be chosen 
With the inner and outer scalar distribution radii, Vg and Rq, getting 



as Tg, Ts and Rq 



^^The gravitational attraction force mrg/2r^, acting on a test particle with mass m inside a lacuna 
due to a central fracture, is superseded eventually by the universal attraction force mvl/r due to halo, 
which should dominate at r > max(i?Oi ^g/w^). 
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Figure 5: Representative irregular solutions X: solid lines - numerical results for interfer- 
ence term Xi/5i, with 5i as in Figure |4j dashed line - piece-wise analytical approximation, 
with So = 0.85 and tq = 2.0. 



closer, the three-component structure of a lacuna becomes less prominent. Nevertheless, 
the asymptotic behaviour and gravitational confinement in a lacuna should still survive 
(at least in a parameter region). This is because the asymptotic properties of the solution 
are determined exclusively by the exponential term in the r.h.s. of the vacuum scalar-field 
equation (77) or (78). 



4 Interpretation and applications 
4.1 Energy content 

Static space structures The preceding results were obtained exclusively in the geom- 
etry framework without any recourse to DM. To the latter end let us reveal the energy 
content of the static space structures found previously. Consider an isolated gravita- 
tionally bound system. Let T^j, be its total bare energy-momentum tensor, with gravity 
included only in the minimal fashion through metric. To incorporate the proper gravity 
contribution define the net energy- momentum pseudo-tensor T^^, = T^^ + f^^, with f^y 
being the pseudo-tensor of the gravity itself. At that, f^y is well-known to be dependent 
both on definition and, possibly, the class of coordinates. The same concerns thus T^^. 
According to Tolman [16], with f^y taken as the Einstein pseudo-tensor in the quasi- 
Galilean coordinates, the net gravitating energy/mass of a static isolated system may be 
expressed entirely through its bare energy-momentum tensor as follows: 

Mg = [fl^£x= [{T^.-TD^gd'x. (98) 



In other terms, the net energy density, with gravity properly incorporated, is p = Tg = 
p + Spi, where p and pi {I = 1, 2, 3) are the bare energy density and pressure, respectively. 
Physically this means that gravity "pumps-in" energy equal to the proper work performed 
against pressure. Having been established in the particular spatial coordinates, Mg is to 
be prolonged to the arbitrary ones as a scalar. For validity of the isolation assumption, 
the ever-present cosmological constant A, which is paramount for the Universe as a whole, 
is to be excluded from p for an isolated system. (Otherwise this would result in a double- 



counting when considering the evolution of the Universe.) In view of (26) we get most 
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generally 



(99) 



with the net gravitating energy determined thus entirely by the i?Q-component of the 
Ricci tensor. 

In the reciprocal gauge, by means of (44) and the first equation (41 ) (with 1/h = a = A 
and a substituted by the spatial scalar In A) we get 

2Rl = {AC{ln A)'y/C = -g^^Vk^i In a. (100) 

Here is a spatial component of the four-dimensional covariant derivative. This relation 
is now valid in the arbitrary spatial gauge/coordinates. This gives 



In view of ( 15 ) (with X substituted by In a) one has 

v^/'VfcV/ Ina = dk{^g^^di\na), 



(101) 



(102) 



and the three-dimensional Gauss theorem reduces (101) to the integral over the remote 
two-dimensional surface S as follows: 



gg'^^dilnadSk, 



(103) 



with ^/—g = y/abc. 

Partite Mg onto the contributions of the non-gravitational matter and the systolons, 
respectively, and M^, including the proper gravity energy: 



Mm + Ms 



= I {(Tm'o - TrZ) + {Tsl - TsD) ^g Sx. 
Put by default for the net partial energy densities: 

7^ —T^ — T"^ n — T — T 
In view of (38) the systolon net partial energy (with Vg neglected) is 

Ms = -vln] J y/^g'^'VkViXd'x, 
and the three-dimensional Gauss theorem gives then 



2 2 



gg'^'diXdSk. 



(104) 



(105) 



(106) 



(107) 



Treated as general covariant scalars the expressions above determine in UBG the 
energy content of a static space structure through V ■ Vina and V ■ VX. While the 
first term is the same as in GR with a (free massless) scalar field, the second term is 
peculiar to UBG. It is produced due to the unimodularity constraint missing in GR^ In 
particular, in the case Ina = v^X (modulo an additive constant) there follows Mg = Mg, 
Mm = signifying a matterless, pure systolon static space structure. 



Physically, the scalar-graviton field may be attributed to (singularity of) modulus as a scalar source. 
This may take place even in vacuum. In GR, to reproduce Mg one should couple the scalar field directly 
to a matter scalar charge. Thus, though the solution in the two theories is formally the same, its physics 
content differs 1241. 
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Dark fractures In what follows we adopt for a dark fracture the truncated real form, 



with a singularity only at r = (see Sec 3.2). Everywhere, but for the origin, there fulfills 



Rq = and V ■ VX = 0. Thus, the spatial integrals are saturated in this case exclusively 
at r = by the point-like singularity. The latter being in fact known only implicitly, 
a spatial integral may be substituted via the three-dimensional Gauss theorem by the 
respective integral over the remote sphere. Having consistently defined the energy in the 
quasi-Galilean coordinates, one can use, just for calculations, the polar coordinates, with 
di = S^dr, y/—gdSk = and g"^"^ = —A. The net gravitating energy of a fracture 

is then as follows: 

M, = AnKli^jrj = r,/(2G) > 0, (108) 
with Tg > {vf > 0) imposed. Likewise, the systolon net contribution is 



= iTiKlvsVf^Jl - i^jrf = Vsrs/{2G), (109) 



or otherwise, Mg = VgTjf^l/v'j — IMg. The net contribution of the non-gravitational 
matter is then 

M^ = Mg-M,= (l- VsVf^l/u} - l)Mg. (110) 

Imposing by default the requirements Mm > and Ms > (r^/ = 1) one can envisage 
two following extreme cases. 

(i) BHs: Vf = i^/max = 1, with M^ = and Mg = M^- Here Vg is determined 
conventionally by the matter net energy, Vg = 2GMm- 

(a) Vacuum dark fractures: uj = lyfmm = '^sl \[^^^rt^-, with = and Mg = Mg. 
Here Vg = VsTs, with Vg ^ at Vs ^ 1. This accords with In a/ = v^Xf (modulo a 
constant) as a generic condition of the absence of matter, 



Dark halos For a dark halo the solution is regular, and according to (38) the systolon 
net energy density is well-defined as 

Ps = -2A, = -2Aoe-^. (Ill) 

Under the attractive effective potential (Aq < 0) this energy density is positive-definite. 



with ps = only at Aq = 0. In view of (76) and the first equation (78) one has 



Substituting X = Xh one then gets in particular the halo energy density profile as follows: 

P^-P^\ l/(3r2) + 0(l/r^/2), T » 1, ^^^^^ 
with Po standing for the central energy density: 

Po = -2Ao = Qvl-§^. (114) 

ito 



The gravitational attraction of the "normal" fractures (^7/ = 1, ^'s > 0) is enhanced, Mg > Mm- 
Admitting rjf = —1, relaxing thus requirement M^ > 0, one would get the "anomalous" fractures 
screening the matter, Mg < Mm- 
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Figure 6: Normalized energy density profiles p{t)/pq: the difference between a profile 
and the reference one Pref/Po = 1/(1 + ''"^)- Bold solid line - numerical result for the 
soft-core profile Ph/po] dashed line - piece-wise analytical approximation; thin line - the 
hard-core profile p/po = l/3r^ for exceptional solution. 



Thus, a dark halo naturally enjoys the soft-core energy density profile, with Rq being 



the halo core radius. At that, ph/po closely reproduces the dark modulus (88). For 



comparison, the exceptional solution X{t) = ln3r^ results in the hard-core profile 



2 



reproducing, accordingly, hard-core modulus (85). Asymptotically, p/t(r) oscillates around 
p with attenuation, approaching p at r ^ 1. The normalized difference, (p — pref)/pQ, 
between a profile and the reference one, pref = Po/(l + t^), is shown at the moderate r 
in Figure [6] 

It is seen, in particular, that the soft-core profile, ph, decays at r < 2 faster than the 
reference one, both having the same behaviour near the origin]^ 
The halo net gravitating energy interior to r is 

At r ^ Rq, one has Mh^{r) = STcv'^K^r = v^r/G, with the respective energy for X rising 
linearly identically. 



To clarity the physics content of a dark halo note that according to (37), (76) and 



(111 ) the exceptional solution in the leading i;<j-approximation results in the bare energy- 
momentum tensor with 

Ipo 
3r2 

In turn, this results in the positive-definite net energy density p = po + ^Pn = Pr = 
Po/(3r^) > 0, being attributed entirely to tensor gravity. As for the regular solution Xh 



Po=Pe=P^ = 0, Pr- = ^^- (117) 



{X'f^\r=o = 0), its bare energy density (37) near the origin is negative, while the net one 



(111) is nevertheless positive-definite everywhere (Aq < 0). The same concerns the bare 



radial pressure (37). At r ^ 1 the regular solution behaves as X, with the net energy 
density decaying like 1/t^, too. Thus, the dark halo is a gravitationally tightly bound 
structure, with the bulk of its net energy provided by tensor gravity to form the halo. 



^''In the conventional DM approach p^e/ is due to the so-cahed pseudo-isothermal sphere, while p 
corresponds to the true isothermal one (see, e.g., [47]). 



29 



Figure 7: Normalized energy density profiles p{t)/pq: solid lines - numerical results 
for interference term pi/{po6i), with 6i as in Figure 4; dashed line - piece-wise analytical 
approximation, with 6q and tq as in Figure 5. 



Dark lacunas For a dark lacuna one has Xi = Xf + Xh + Xi = Xf + Xheff- The 
contribution to the net gravitating energy due to Xf is accounted for by the surface 
integral. For energy density of the effective dark halo, Pheff = Ph + Pi, Eq- (112) gives ph 
as before and the interference contribution as follows: 

n/n /l/r-4r + 0(r3), r < 1, 

Thus, Pheff comprises some cuspy 1/r-correction compared to the soft-core profile ph- For 
the normal fractures {rif = 1) the interference is constructive. Si > 0, with the effective 
halo net energy near the origin increasing. The numerical results for the normalized 
interference term, pi/{pQSi), are shown in Figure [7j The interference contribution to the 
lacuna net energy interior to r is 

f \ TO 2 /r2-2r4 + 0(r6), r < 1, 
Mi^{r) = l2nVsngrA^^^^^ (119) 

The above results are obtained in the perturbative fashion. For the realistic lacunas 
the picture may become more complicated in detail though its salient features should, 
conceivably, survive. 



4.2 Rotation curves and equivalent DM 

Rotation velocity Let = dx^/dr, U ■ U = g^yU^U" = 1, be the four- velocity of 
a test particle with respect to its proper time r. If the particle interacts gravitationally 
only in a minimal fashion through metric, its If^ for the free motion satisfies the geodesic 
equation, dU'^/dr + V^JJ^U^ = 0, with F^^^ standing for the Christoffel connection. 
For the circular rotation in a static spherically symmetric metric one has in the polar 
coordinates U"^ = 0. Besides, one can put 6 = 7r/2, with = dU^/dr = 0. The 
remaining non-zero components are Lf^ = dip/dr and = dt/dr. The ratio u = 
W^/U^ = dif/dt is thus nothing but the angular velocity with respect to the observer's 
time. The visible circular rotation velocity is as follows: 

u = ^/a^d^/dt = ^/a^U'^/U^. (120) 
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For the spherically symmetric metric (36) one has Fqq = a' / (26), r^^|e=7r/2 = —{cr )'/ (26) 
and = 0, so that the equation of motion gives dV/dr = -(rSo(t/°)2 + r;^(?7^)2) = 0, 
and thus 

(f/^)V(f/°)2 = a'/icr^y. (121) 
Altogether, the visible rotation velocity squared is 

= a' / \n{cr^y . (122) 

Otherwise, the rotation velocity with respect to the particle proper time looks like 

V = udt/dT = uU° = Vcr^U'^. (123) 



Accounting for f/ ■ t/ = a{U°f - cr'^{U'^f = 1 and (121 ), one then gets V = vj^l - ' 



where 



(lna)'/(lncr 



2\i 



(124) 



In the above, v = uj^fa = V/\/l + is the visible rotation velocity accounting for 
the gravitational deceleration of time. For consistency, v < 1. Being gauge invariant, 
all these expressions are valid in the arbitrary radial coordinates. The choice of the 
expressions for the rotation velocity remains though convention dependent. By default, 
we choose v. Particularly, in the astronomic coordinates (c = 1) one gets conventionally 
t;2 = r(lna)'/2. In the non-relativistic weak-field limit, we are interested in, the difference 
between the definitions becomes irrelevant 
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Rotation curves In the leading t;5-approximation one can put c = 1, so that 

= ra'/2 = {T/2)da/dT. (125) 

The total velocity squared in a dark lacuna is thus the sum of three components, vf = v'j + 
v\ + vf = vj + vl^fj, where the fracture, halo and interference contributions, respectively, 
follows: 




|r4 + ir6 + 0(r8), r < 1, 



l + 0(l/v/^), r»l. 



with vf being generally sign-indefinite. It is seen that the RC profile Vheffir) for the 
effective halo gets flat only asymptotically. At that, the exceptional solution X would 
result in the precisely fiat profile v = Vhoo = Vs, around which all the profiles Vheffij) 
oscillates with attenuation, approaching at r ^ 1. The respective results are shown in 
Figures |8 10 Evidently, the net velocity profile is in reality far from being exactly fiat. 



^^In passing, choosing u one would get that for BHs (a = 1 — r^/r, c = 1) in all the aforementioned 



radial coordinates (see Sec. 3.2), but for isotropic ones, there exactly fulfills the third Kepler's law: 
~ r^, with T being the rotation period. 
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Figure 8: Normalized RC profile v{t): bold solid line - numerical result for v][/v1 
practically coinciding with analytical approximation (dashed line); thin line - {v'^ /v1 = 1) 
for exceptional solution. 



0.02 - 




Figure 9: The asymptotic of the normalized RC profile v{t): solid line - numerical 
result for {v\/v1 — 1) practically coinciding with analytical approximation (dashed line). 




Figure 10: Normalized RC profile v{t): solid lines - numerical results for interference 
term vf/{v'^6i), with 6i as in Figure 4; dashed line - piece- wise analytical approximation, 
with So and Tq as in Figure 5. 
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Equivalent non-relativistic DM Let us interpret tlie previously found energy den- 
sity and rotation velocity profiles in terms of the non-relativistic DM, the latter being 
conventionally used in confronting with astronomic observations. For the non-relativistic 
DM halo in the fiat space-time (a = c = 1) the Newton's dynamics results in the rotation 
velocity Vd determined through 



vyr = GMd<ir)y, 

where Md<{r) = An Jq Pd{r)r'^dr is the DM energy interior to r, with 
energy density. The latter should thus satisfy 



{rvjy 



Pd 



AttG r2 



(127) 
being the DM 



128) 



Imposing vj = vl^ff = 
equivalent DM as follows: 



+ fj and accounting for (125) we get the energy density for the 



Rq dr \ dr 



(129) 



where aheff = o.h + on- Further, because in the second line of (78) the term 1/r due 
to fracture drops off, one gets hereof aheff = v^Xheff, and (112) for Xheff gives pd = 
Pheff = Ph + Pi- The profiles Pd^T) oscillate with attenuation around p approaching the 
latter at r ^ 1 (cf., e.g.. Figure 9). 

Thus in the context of RCs, the previously found energy content of the dark halo co- 
incides in a consistent manner with the non-relativistic DM interpretation. At that, RCs 
admit a complementary approach: either directly in terms of the static space structures 
in the geometry framework, without any recourse to DM, or in terms of the equivalent 
non-relativistic DM in the framework of the Newtonian dynamics. This justifies treating 
the coherent systolon field as DM. 

For more clarity let us present pd as follows: 



Pd = (Po/6)x d Xheff I dx 



where x = Po is given by (114) and there roughly fulfils Xheff — ^h, with 



Xh 



1 HVx') 



so that 



Pd ^ Po 



1, 

XV3, 



X>1, 
X< 1- 



(130) 



(131) 



(132) 



(Note in parentheses that x is nothing but the radial harmonic coordinate for the dark 



halo, with the approximate scalar-field equation (78) being d^Xh/dx^ 
Also, under the choice of the unit of length Iq 
(53) smoothly matches at rj 



< r < 



Rq the fracture harmonic coordinate 
i?o with X for halo.) Evidently, it is only the 
relatively rapidly varying with x part of the scalar-graviton field in halo which matters 
for the equivalent DM. At that, varying slower than logarithm (though, possibly, larger) 
part of the field may result in the equivalent DE through the potential Vs{X) (omitted 
here). The spherical symmetry being not crucial, we expect that this is a generic property 
of the systolon DM and DE. 
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We believe that such a coherent scalar-graviton field is of principle importance only 
within the gravitationally bound systems, such as galaxies and clusters of galaxies. Be- 
yond them the specific nature of the interior DM may become less important, with the 
Universe as a whole being described effectively by the conventional ACDM model (or 
a variation of it). This would result in the apparent two-component DM and could, 
perhaps, help unravelling the so-called core/cusp DM problem |47j , 



Galaxy DM: coherent scalar field vs. continuous medium To refine the nature of 
the equivalent DM let us consider an isotropic continuous medium with the conventional 
energy-momentum tensor 



flU 



(133) 



where pd and pd = Pd{Pd) are, respectively, the proper energy density and pressure of the 



putative medium, and f/^ {Ud ■ Ud = I) is its local four- velocity. According to (26), in 



the comoving frame, Ud = (1,0,0,0), (assumed to coincide with that of galaxy) in the 



metric (36) there fulfills: 



^0 = A(aPrf + (« + 2)Pd 



2k 

Imposing in the weak- field (|a — 1| ^ 1) approximation 

Pd + 3pd = Pheff, 



(134) 



(135) 



i.e., equating the net energy densities one reproduces the same Rq for the medium and the 



effective halo. With largely determining In a through (100), and In a, in turn, largely 



describing RCs through (124), this ensures the equivalence between the two descriptions. 



the coherent-field and continuous- medium ones, in the non-relativistic RC context. Under 
(135) such an equivalence fulfills irrespective of the equation of state pd = Pd{Pd)- In 
particular, at pd ^ Pd one recovers the condition of the preceding paragraph, pd = 
Pheff, for the cold/warm DM. As far as being determined largely by the Newtonian 
potential. In a, the same statement concerns other gravitational manifestations in the 
non-relativistic weak-field approximation. 



Soft-core DM halos The observational data on the DM-dominated (late-type LSB 
disk and gas-rich dwarf) galaxies are consistent with the cored energy density profile for 
the DM halos as follows: 

where pc and Rc are two free core parameters (see, e.g., [H]). The presence of a mild 
cusp is still admitted by the data. Such a profile leads to the asymptotic constant 
rotation velocity v^a = {A-kG pcR^Y^"^ ■ This behaviour is in qualitative agreement with 
that for the previously found one-parameter dark halo profile ph{r/RQ), which results in 
the asymptotic constant velocity Vhoo = (47rGpo-Ro/3)^^^- But now the parameters po and 
Ro are, in fact, not independent ensuring the fixed universal Vhoo = Vs- To disentangle the 
parameters one should use a three-parameter lacuna profile accounting for the fracture- 
halo interference. The latter gets significant with the interference parameter 5i = 0{1), 



29 



Besides, as an incoherent component of DM built of systolons there might serve the (near) vacuum 



mini-fractures (see Sec. 4.1) distributed all over the Universe. 
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or with the scalar radii for fracture and halo related as ~ VgRo, Vg 1, what is not 
quite unrealistic. Other factors such as lacuna asphericity or rotation, influence of the 
scalar-field potential, influence of the distributed non-gravitational matter, etc, may also 
be of importance, ^^1^ 

With these caveats, the dark lacunas consisting of a stabilizing super-massive dark 
fracture at the origin surrounded by a dark halo might serve as a prototype model for the 
galaxy DM frames, to be supplemented ultimately by the distributed non-gravitational 
matter. With Kg = Mp/V^ = 2.4 x 10^^ GeV, Mp = 1/y/G standing for the Planck 
mass, the asymptotic rotation velocity in galaxies foo ~ 10^'^ (in units of the speed of 
light) at its face value results in Vg = y/2Ks/Kg ~ 10"'^. This implies that the mass scale 
appropriate to the scalar mode, Kg ~ 10^^ GeV, is to be of the order of GUT scale. Could 
it be more than just a coincidence, with a common origin of the two scales (if any)? 



Long-distance gravity modification: vacuum vs. Lagrangian In the end let us 

present several comments concerning conceivable long-distance gravity modification in 
the context of galaxy DM. 

(i) Spontaneous vacuum modification The present approach to galaxy DM com- 
prises an explicit modification of the gravity Lagrangian at a fixed ultraviolate mass 
scale Kg ~ 10~^K,g. The effective theory is to be valid up to the high scales fi < Kg. An 
infrared parameter Aq (or, equivalently, a long-distance scale Rq) appears at random due 
to spontaneous breaking of a global symmetry within a lacuna. The randomness of the 
infrared parameter insures more flexibility of the theory in the context of galaxy DM. 

(a) Explicit kinetic Lagrangian modification The attempts at the explicit infrared 
Lagrangian modifications in the context of galaxy DM are numerous. First of all, one 
may mention the so-called /(i?)-gravity, without or with a proper modification of mat- 
ter Lagrangian (cf., e.g, [IHl ED])- In this case, having no specific degree(s) of freedom, 
DM is just mimicked by the modification of the long-distance tensor gravity]^ A re- 
lated approach is given by the scalar-field theories with a non-canonical kinetic term, 
supplemented ultimately by a repulsive potential in Lagrangian (cf., e.g., [511 [52]), etc. 



Moreover, one may envisage the case wth a ghost quadratic kinetic term,^^ By construc- 
tion, so modified Lagrangians are given by some functions /{R/fij) or /(Vcr • Va/fij), 
etc., depending explicitly on a fixed infrared mass scale fii -C Kg. Ultimately, this would 
imply that a more fundamental theory should settle down already on the relatively low 
scales fi > fip 
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•^"In particular, despite the fact that LSB galaxies are DM-dominated, the account for their star disks 
may be important for RCs near the origin. 

•^^From the DM point of view the dark lacunas may rather be considered as scalar "lumps" in asymp- 
totically non-flat metric. This is to be contrasted with alternative attempts at treating galaxies in 
the framework of GR with a scalar field by means of scalar lumps in asymptotically fiat metrics (cf., 

e.g., m)- 

•^^ At that, a looked-for putative matter is rather "missing" than the dark one, being beyond the reach 
of direct searches. 

■^^For the scalar-field ghost condensation as an alternative to DM in the context of the Universe cf, 
e.g., [53]. 

•^^It goes without saying that UBG could a priori admit, if desired, an arbitrary infrared modification 
of the kinetic Lagrangians both for the tensor and scalar modes. Presently such modifications are left 
aside. 
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5 Conclusion and prospects 



The Unimodular Bimodal Gravity (UBG) is a theoretically viable development of UR and 
GR with a scalar field. It retains the principle ingredients of GR - the general covariance 
and masslessness of the transverse-tensor graviton. Beyond GR, the theory comprises in 
metric a propagating compression mode - the scalar graviton/systolon. With the latter 
treated as DM, UBG presents a unified description of the (tensor) gravity and DM. The 
appearance of a physically well-motivated scalar field is the crucial point of the theory. 
In its reduced version, with the spontaneously broken dynamical global symmetry, the 
theory, being rather restrictive, is apt to result in a number of definite predictions. In 
particular, in the static spherically symmetric case it predicts peculiar space structures - 
the dark lacunas - consisting of a compact singular dark fracture (a scalar-dressed BH) 
at the origin surrounded by an extended soft-core dark halo. Enjoying the property of 
gravitational confinement, with the logarithmic potential of gravitational attraction at 
the periphery, the dark lacunas ensure asymptotic fiattening of RGs and may serve in 
cosmology as the DM frames for galaxies. 

The scalar graviton/systolon physics presents conceivably a perspective field of the 
future investigations, both theoretical and phenomenological. In particular, the infiuence 
of the (near) massless scalar field may drastically change the structure of the GR BHs 
both at their event horizon and at asymptotic due to appearance of dark lacunas. Further 
studying the latter ones, as well as dark fractures and halos, including their more subtle 
aspects such as a putative asphericity or rotation, influence of the scalar-graviton poten- 
tial and the distributed non-gravitational matter, application to various types of galaxies 
as well as to galaxy clusters is in order. The application of the theory to evolution of 
the Universe is likewise urgent to verify the theory (if any) and unravel ultimately the 
mystery of DM and DE. 

Acknowledgement Thanks are due to I. Yu. Polev for assistance with numerical cal- 
culations. 
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